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ABSTRACT 

The current sample of high-redshift Supernova Type la, which combines resuhs from 
two teams, High-z Supernova Search Team and Supernova Cosmology Project, is ana- 
lyzed for the effects of weak lensing. After correcting SNe magnitudes for cosmological 
distances, assuming recently published, homogeneous distance and error estimates, we 
find that brighter SNe are preferentially found behind regions (5-15 arcmin radius) 
which are overdense in foreground, z ~ 0.1 galaxies. This is consistent with the inter- 
pretation that SNe fluxes are magnified by foreground galaxy excess, and demagnified 
by foreground galaxy deficit, compared to a smooth Universe case. The difference 
between most magnified and most demagnified SNe is about 0.3-0.4 mag. The effect 
is significant at > 99% level. Simple modeling reveals that the slope of the relation 
between SN magnitude and foreground galaxy density depends on the amount and 
distribution of matter along the line of sight to the sources, but does not depend on 
the specifics of the galaxy biasing scheme. 

Key words: gravitational lensing - supernovae: general 



1 INTRODUCTION 

The effects of weak gravitational lensing by the large-scale 
structure have been detected in several samples of high 
redshift QSOs, intermediate redshift galaxies, and BATSE 
GRBs. In the case of point sources, QSOs and GRBs, 
weak lensing manifests itself as angular (anti-)correlations 
between these sou rces and foreground in h omogeneously 
distributed mass ("Willia ms fc FrevI l2003l: iB enitez et al.' 
boot; Wilhams & Irwin 19981: iBenitez fc Martine z-Gonzalez 
|l99S), while in the case of galaxies we ak lensing i s dete cted 
through its coherent shear effect (see iRefreeien (|200^ for 
a recent review). In principle, there is another, more direct 
way of detecting weak lensing, which uses fluxes of standard 
candles. If the observed magnitudes of standard candles are 
corrected for cosmological distances then the effect of lensing 
can be seen: brighter sources will lie behind regions of mass 
density excess, while fainter ones will have mass deficits in 
their foregrounds. 

The best example of cosmological standard candle, Su- 
pernovae Type la (SNIa) have been extensively observed 
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extinction have been studied theoretically and observation- 
ally, and have either been corrected for or shown to be 
small. Weak lensing, another nuisance e ffect has been ad- 
dress e d theoretically by several authors (lAman ullah ct alj 
120031 iMetcalj Il999l: iHolzl Il998l : IWambseansset al. .199% 
and found to be u nimportant given the cur rent uncertain- 
ties. For example, 'Wa mbseanss etlil] 1^93) used ray trac- 
ing through cosmological simulations and found that the 
lensing induced dispersions on truly standard candles are 
0.04 and 0.02 mag at redshift z = 1 and z = 0.5, respec- 
tively, in a COBE-normalized cold dark matter universe with 
Qm = 0.4, Qa = 0.6, Ho = 65km/s/Mpc and ag = 0.79. 
These are small variations compared to the current errors 
which are > 0.2 mag. Even though weak lensing effects are 
estimated to be small for Zs < 1, they are predicted to be 
non-negligible for higher redshift sources, so it is not surpris- 
ing that the highest redshift SNIa, SN199 7ff at Zs =_l-755 
has been examined by several authors ( | Benitez et al.l l2002l : 
iMoertsell et~alll200lt iLewis fc Ibatall200ll) for the effects of 
weak lensing due to galaxies along the line of sight. 

Present day high-z SNIa samples are dominated by 
lower redshift SNe, and so have not been examined for the ef- 
fects of lensing. The main goal of this work is to determine if 
the observed fluxes of the cosmologically distant SNIa have 
suffered significantly from lensing induced (de-) amplifica- 
tions. 
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2 DATA 

The largest homogen eous compila t ion of SNIa has been re- 
cently published by iTonrv et all 1)200.^^ : Table 15 of that 
paper contains 74 SNe at Zs ^ 0.35 . The authors use 
four different light curve fitting methods (MLCS, Ami5(S), 
modified dml5, and Bayesian Adapted Template Match) to 
estimate distances to SNe. The final quoted distance is the 
median of the estimates of the four individual methods, and 
the uncertainty is the median of the err or of th e contribut- 
ing methods. The analysis presented in iTonrv e t al. (2 0o3) 
yields values of the global cosmological parameters; if a fiat 
model is assumed, then Qm = 0.28 and JIa = 0.72. We use 
these values in all the analysis of the present paper. 

As t racers of foreground mass density we use APM 
galaxies dlrwin et al. 1994'). APM provides near full cover- 
age of the sky in the Northern and Southern Hemispheres, at 
|6| > 20°. In our analysis we use only the central r = 2.7° of 
APM plates. Since the plate centres are separated by 5°, 
there exist small portions of the sky that are not covered by 
any plate. As a result of these cuts, only 55 of the 74 SNe 
lie on the usable parts of APM plates. 

The median redshift of the 55 SNe is 0.47.^ Since most of 
the SNe have rather low redshifts, care must be taken to en- 
sure that galaxies are foreground to the SNe. Furthermore, 
because SNe span a large range of nearby redshifts, from 
Za = 0.35 to 1.2, the optimal lens redshift zi will depend on 
Zs much more compared to a typical case where sources (say, 
QSOs) are at Zj, ~ 1 — 3 and so the redshift of optimal lenses 
is roughly independent of Zs- In our analysis we adjust zi for 
each SN source by selecting the appropriate limiting appar- 
ent magnitude , mag ggi um for APM galaxies on red plates. 
iMaddox et al.l il99(fl gives an empirical expression for the 
median redshift Zmed of a galaxy sample with a given faint 
magnitude flux cutoff. This median redshift can be equated 
with the optimal lens redshift zi, and hence the magnitude 
limit of the foreground galaxies can be determined for every 
SN separately. However, there is a small catch. For Zs = 0.4 
optimal zi = 0.174. The galaxy redshift distribution whose 
median redshift Zmed ~ 0.174 has a considerable tail extend- 
ing beyond z — 0.4. To avoid the problem of source/lens 
redshift overlap we use Zmed = zi/2, where factor of 2 was 
chosen arbitrarily. We explore the dependence of the results 
on this factor in Section |1| 



3 ANALYSIS 

Around every SN we draw a circle of radius 9 = 10', and 
count the number of galaxies, ngai,D, in the appropriate 
magnitude range. This number is compared to the average 
number density in control circles, {ngai,R}- Fractional galaxy 
excess is Sugai = ngai,D/{ngai,R) — 1. Control circles are con- 
fined to the same APM plate as the SN, and to the same 
distance from the plate centre as the SN (to reduce the ef- 
fects of vignetting); however, scattering the control circles 
randomly on the plate does not change the results signifi- 
cantly. For each SN we also calculate Ny/N, where A''> is 
the number of control circles, out of total A^, that have less 

^ SN1997ff at Zs = 1.755 is not in our sample: it fell in the cracks 
between the APM plates. 
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Figure 1. Absolute magnitudes of SNIa (modulo a constant off- 
set) versus the foreground galaxy density. Magnitudes were ob- 
tained assuming f2„i = 0.28, and Qa = 0.7 2, and extinction, K- 
corrected distances from lTonrv et aU (2003). A'> /N is a measure 
of the number density of foreground galaxies in circles of radius 
9 = 10' around SNe (see Section|21for details). There are a total of 
55 sources (with magnitude errors), but only 50 (filled points) are 
used in the analysis. The 50 points are grouped into three bins, 
whose horizontal size is shown as thick horizontal "error-bars". 
The corresponding vertical error-bars show the deviation of the 
mean of the points in each bin (the rms is 4 times larger), and the 
intersection of the thick lines are the averages of the Magsjv of 
SNe in each bin. The thin slanting line is the best-fit to 50 SNe, 
and has a slope /3 = —0.373. The significance of the correlation is 
> 99%. 

galaxies in them than the circle around the SN. In other 
words, A'^> /N is the rank of the SN circle among its con- 
trol 'peers'. If SNe are randomly distributed with respect 
to the foreground galaxies, then average (A''> /N) = 0.5. If 
SNe have an excess (deflcit) of galaxies in front of them then 
their {A*'> /N) will be greater (less) than 0.5. Analogous to 
the medians being more stable than averages, Ny/N rank 
statistic is more stable than Sngai- 

Figure shows absolute magnitudes, Msn, and Ny/N 
ranks of 55 SNe found on APM plates. The effect of fiux 
dimming due to cosmological distances has been taken out, 
i.e. all the SNe have been 'brought' to the same redshift; the 
constant magnitude offset on the vertical axis is irrelevant 
for this work. There are two SNe whose magnitudes make 
them > 3ct outliers, SN19970, and SN1997bd represented by 
empty circles with dotted line error-bars. We exclude these 
from our analysis.^ 

2 Of our 55 SNe, SN1997b d has the h i ghest host-galaxy ex- 
tinction, (A„) as quoted in iTonrv et al.l i2003^ . SN19970 and 
SN19 94H were excluded from the primary fit (Fit C) by the anal- 
ysis of lPerlmutter et al.ltT999l) . We do not exclude SN1994H from 
our analysis, but if we did it would improve the trend, as its co- 
ordinates on Fig.0are (0.427; -0.606). 
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Figure 2. Histogram of SNIa magnitude errors , SMcg^ = 
5 (51og (d). Distance errors, Slog(d) were derived bv lTonrv et alj 
J2003r . and listed in column 9 of their Table 15. The solid his- 
togram is for 55 SN; the dashed line represents 31 SN from 
the High-2 Supernova Search Team. The 3 outlier SNe with 
SMsN > 0.45 are not used in the present analysis. 



The dis t ance (or, magnitude) errors, as estimated by 
iTonrv et all ^2003^ for the 55 SNe are shown in Fig. ^ 
as the solid line histogram. The dashed line represents the 
subsamplc of 31 SNe from High-z Supernova Search Team 
(Tonrv e t al. 20Q3i), whose errors appear to be generally 
smaller than those of the Supernova Cosmology Project 
ijPerlmutter et a^]ll999^ . We use SNe from both the groups, 
but exclude three outliers in Fig.|51 whose SMsn > 0.45. In 
Fig. these are represented by empty circles with dashed 
line error-bars. Thus, we exclude a total of 5 SNe from our 
analysis, leaving us with 50, shown as the solid points in 
Fig.El 

These 50 SNe exhibit a relation between their Msn 
and A''> /N, in the sense that brighter SNe have an excess of 
galaxies in their foregrounds. For illustration purpose only, 
we bin the 50 SNe into three bins; Fig. shows the extent 
of the bins and the deviation of the mean of the SNe mag- 
nitudes in each bin as thick lines. 

The best-fit line to the 50 SNe in Fig. has a slope 
/3 = —0.373, and is shown as a thin slanting line. This fit 
does not include magnitude errors. To include the errors we 
do the following. We calculate the best-fit slope for 10,000 re- 
alizations of the data, where each data point's SN magnitude 
is replaced by a randomly picked magnitude from a Gaus- 
sian distribution centred on the actual magnitude value, and 
having width equal to the quoted error. This procedure cor- 
rectly incorporates the information contained in the errors, 
and produces a distribution of best-fit slopes, which is shown 
as a solid line histogram in Fig. 01 This distribut ion shows 
that the median best-fit slope is — —0.372, while (3 = 
(i.e. a case of no correlation) is ruled out at 99.80% confi- 
dence level. Had we used 53 SNe (i.e. had we not exluded 



1 — ' — ' — ' — I — ' — ' — ' — I — ' — • — ' — I — ' — ' — • — \ — • — ' — ' — I — r 




-0.8 -0.6 -0.4 -0.2 0.2 

|8 

Figure 3. How do the SN distance errors affect the best-fit slope, 
/3? The histograms show the distribution of derived (5 values in 
10,000 realizations of the data. The solid (dashed) line histogram 
represents results using 50 (53) SNe. The short vertical lines at 
the bottom of the plot indicate the medians of the distributions. 

the 3 SNe with large distance errors), the median best-fit 
slope would gave been f3 = —0.336, while /3 = would have 
been ruled out at 99.37% confidence level. The correspond- 
ing distribution of P slopes is shown as the dashed line in 
Fig.H 

Next, we determine the likelihood of this relation aris- 
ing by chance. To that end, we estimate the significance of 
the relation in two ways. First, we assign random positions 
to SNe (keeping these control SNe on the same APM plate 
as the original source), and redo all the analysis. We repeat 
this 1000 times, and in 8 cases we find (3 ^ —0.373, im- 
plying statistical significance of 99.2%. Second, we take the 
list of observed Msn values and randomly reassign them 
to observed SN sky positions. 10,000 randomized SNe sam- 
ples are created in this fashion, and only 0.5% of these have 
fH ^ —0.373. Based on these two tests we conclude that the 
significance of the MsN-Ny /N relation is better than 99%. 

These results are consistent with weak gravitational 
lensing, which would amplify SNe found behind more nearby 
mass concentrations, as traced by APM galaxies. Alterna- 
tively, the results could be due to the action of Galactic 
dust, which will obscure certain directions of the sky making 
galaxies less numerous and SNe fainter. We consider Galac- 
tic dust further in Section|3 in Sections|^and Elwe proceed 
on the assumption that weak lensing is responsible to the 
MsN-Ny/N relation. 



4 FURTHER TESTS 

The distribution of SNe points in Fig. Q depends on specific 
choices that we made for certain parameters, in particular 
we chose circles of radius 6 = 10' and galaxy magnitude 
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Figure 4. Control tests: varying the aperture radius, d of circles 
around SNe. The vertical axis shows /3, the slope of the best-fit 
line to the relation of the type shown in Fig The filled cir- 
cles represent cases where galaxies were counted around SN; stars 
symbols represent cases where Galactic stars were used instead 
(i.e. control experiment). The dotted lines denote, approximately, 
the levels of statistical confidence. The dashed vertical line indi- 
cates radius d used in the analysis of Section ISl 



limit such that Zmed = zi/x, where x = 2 (see Section I^J. 
How would the results change if different choices were made? 
In other words, how robust is our result, would it disappear 
had we picked a different set of parameters? 

Figure|l]shows the effect of changing 9. The vertical axis 
is /3, the best-fit to the MsN-Ny /N relation in each case. 
Filled points represent cases where galaxies were counted 
and used to determine the A''> /N rank, while star symbols 
represent cases where Galactic stars were used instead. As 
expected, Suatara do not correlate with SNe magnitudes, and 
the values of the best-fit slopes are near zero. However, be- 
cause the APM star-galaxy classifier is not perfect, some 
'stars' are actually galaxies, which accounts for some signal 
being seen when using Sriatara ■ Horizontal dotted lines mark 
the approximate location of 90%, 95%, and 99% confidence 
levels. These are only approximate because every point in 
the plot will have its own significance level, but because the 
number of SNe contributing to each point is the same in 
each case, and the total dispersion in SNe magnitudes is the 
same, same /3 values have about the same significance, re- 
gardless of 6. We note that for very small 6 the galaxy num- 
bers become very small, and Poisson noise drowns out any 
MsN-Ny /N correlation that might exist, so the upturn in 
the values of the best-fit slope at small 6 is probably not real. 
The dashed vertical line marks the 9 value used in Section |21 
We conclude that significant Msn-N^ / N anti-correlations 
occur only with galaxies and not with Galactic stars (which 
serve as a control sample), and only for 5' ^ 9 < 15'. 

Figure 1^ shows the effect of changing the median red- 
shift, Zmed ~ zi/x of the APM galaxies, or equivalently. 



Figure 5. Control tests: varying the apparent magnitude limit 
of galaxies and the corresponding median redshift Zmed- Similar 
to Fig.|^ see Section]^ for details. 

mag gai lim- The hidden variable which is varied along the 
horizontal axis is x. Since each SN has a different value of 
zi and hence msLg gai lim, depending on its Za, there is no 
unique way of labeling the horizontal axis by using galaxy 
magnitudes, or redshifts. We label that axis by assuming 
Za = 0.47, the median of the SNe redshift distribution. At 
the top of the plot we show how galaxy magnitude limit on 
the horizontal axis translates into the median galaxy redshift 
(for a specific case of Za = 0.47). We see that significant anti- 
correlation between Msn and A''> /N of foreground galaxies 
occurs for galaxies with Zmed Ss 0.12. Using Galactic stars 
instead of galaxies produces no significant results. The value 
we used in Sectionj^is shown with a vertical dashed line, and 
corresponds to z^ed ~ 0.1. As in Fig. 2]the upturn in the 
best-fit slope values at bright mag gai Um is probably due to 
Poisson noise. The range of redshifts of APM galaxies that 
act as the best lenses for SNe are ~ 0.05 — 0.2; more distant 
galaxies show no signal, as expected, if lensing is the correct 
interpretation of the data, because more distant galaxies are 
either too close to SNe in redshift or are actually at the same 
z. 

A number of other tests have been carried out as well. 
For example, instead of using a Zs-dependent mag gai iim, 
we tried fixed values of 17.5, 18.0, 18.5, which gave, (5 — 
—0.35, —0.31, —0.27, respectively, comparable to, but some- 
what smaller than those seen in Fig. El This is not surprising: 
Zs-dependent galaxy magnitude limits pick optimal lens red- 
shifts for each source, thus maximizing the observed lensing 
signature. 

We also reran the analysis with subsamples of the entire 
50-source sample. We split the SNe according to the teams: 
20 Supernova Cosmology Project SNIa gave a Msn-N^ /N 
anti-correlation significant at 95.8%, while the correspond- 
ing significance level for the 30 SNe from the High-z Su- 
pernova Search Team is 96%. Because the quality of UKST 
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APM (Southern Hemisphere) plates is higher than POSS 
APM (Northern Hemisphere) plates, and because there is 
some overlap between plates along the equator, we used 
UKST plates whenever possible. Redoing the analysis using 
only the 36 UKST SNe we get the best fit slope (3 = -0.46 
at a significance level of 99.6%, while the 14 POSS SNe had 
P = —0.11 and the correlation was not significant, which 
is not surprising given the size of this subsample. Splitting 
the whole sample into low and high redshift groups we get 
the following results: 25 SNe with Zs < 0.47 have a slope 
(3 = -0.37 at 96.8%, while 25 SNe with with Zs > 0.47 have 
a slope P — —0.23 at 88.9%. These tests suggest that the 
MsN-Ny. /N relation has a physical origin (weak lensing or 
Galactic dust), and is not an artifact arising form one subset 
of the data. 



5 MODELING THE OBSERVATIONS 

In this Section we adopt the weak lensing interpretation of 
the MsN-Ny /N relation. Our simple lensing model does not 
use the individual values of Zs, and the optimal redshift dis- 
tribution of the corresponding lenses. In lieu of these param- 
eters we use the lensing optical depth of the matter traced 
by the galaxies, kq- For a fixed global geometry, kq depends 
on Zb and the redshift distribution of the mass traced by the 
galaxies. We assume that the APM galaxies faithfully trace 
the mass up to some redshift, zt, then 



KO = 0.m Peril (c/ Hq) 



dz. 



(1) 



\Q.m{l + z)'J'+^AY/^Y.crrt{z,Z,) 

In that case, projected fractional mass excess is Srimass = 
ac/ko, where k is the convergence, with respect to a smooth 
Universe (filled beam) for the corresponding line of sight. 
In the weak lensing regime source amplification A ~ {1 -\- 
2k); a = 1 for an unlensed source. The smallest value that 
K can attain is minus the total optical depth, —kq, which 
corresponds to emptying out all the mass along the line of 
sight from the observer to zt- 

In addition to the value for kq, our model has three in- 
gredients, (1) mass distribution, i.e. a probability distribu- 
tion (pdf) for Srimaas, which is related to k pdf, (2) a bias- 
ing scheme, which relates Srimasa to the projected fractional 
galaxy number density excess, Sugai, and (3) the dispersion 
in the intrinsic magnitudes of SNIa's. Given this informa- 
tion we generate synthetic SNIa samples, with 50 sources 
each, then compute Ny/N each SN, and the observed SN 
magnitude, Msn- The slope of the best-fit line to the Msn- 
Ny /N relation, f3 is then used to test how well a given model 
reproduces the observations. 

Because our model assumes the same zi distribution for 
all sources, the value of /3 to compare our models' predic- 
tions to should be the one obtained by using a constant 
mag gal lim ~ 18, i.e. /3 ~ —0.3 (see fourth paragraph in 
Section 

The specifics of the three model ingredients, and the 
associated model parameters are described later in this Sec- 
tion. The goal of the modeling is to determine what set of pa- 
rameters can reproduce the observations, i.e. have /3 ^ —0.3. 
We are particularly interested in what effect biasing has on 
the results. It is often suggested that the large amplitude of 



QSO-galaxy correlations mentioned in the Introduction is, 
at least in part, due to the fact that biasing is not a sim- 
ple linear, one-to-one mapping from Sumaas to Sugai- QSO- 
galaxy correlation function (loqg), galaxy autocorrelation 
(ijJgg) and the matter fluctuation power spectrum (Pk) are 
related by 



LJQG OC b ^LUGG OC bPk, 



(2) 



where 6 is the biasing parameter. This means that uoqg 
probes a combination of 6 and Pu- Thus, for a given Pk, 
i^QG can be enhanced with positive biasing, especially if it 
is non-linear on relevant spatial scales. 

The important difference between weak lensing induced 
QSO-galaxy correlations and the effect on standard candles 
we have studied here is that in the latter case biasing plays 
a minor role. The slope of the MsN-Ny /N relation is most 
sensitive to the mass distribution. Our modeling demon- 
strates this later; first, we described the specifics of the three 
model ingredients. 

5.1 Mass distribution 

5.1.1 The standard mass distribution 

Based on the res ults from cosmological N-body simulations 
J Jain et alJl2000D the shape of the probability distribution 
function (pdf) of k is very roughly Gaussian, but asymmet- 
ric, with the most probable k, which we call Hm being less 
than 0, i.e. most sources are deamplified. The tail of the 
K distribution extends to high positive val ues of k. We use 
these published results of Ijain et alj i2000l) to construct an 
approximate shape for our k pdf: 



p(k) = 



exp ( — [k — K,n]^/2(TKi ), for k > k„ 
exp ( — [/t — /2(7k2), for K ^ hi„ 



(3) 



The pdf is a combination of two half-Gaussians, with two 
widths, CTki and ctkj, which describe the high k and the low 
K sides of the pdf respectively. The location of the peak of 
the pdf. Km is adjusted such that the average k is 0, which 
implies 



(4) 



Km is always negative, since the skewness of the pdf dic- 
ta tes that a^. , > g^ o . From the numerically computed pdfs 
of I Jain et alJ i2000h we estimate that for a flat ilm = 0.3 
model, where the sources are at Zs ~ 0.5 — 1, and the smooth- 
ing scale is ~ 5' — 10', the following approximations apply: 
o'ki/o'k2 ~ 2.6, and ctk^ ~ ko/3.6. Thus we now have a 
realistic k pdf, appropriate for a standard cosmology with 
standard mass distribution. Figure|7[d) shows such a pdf as 
a short-dash line. 



5.1.2 Non-standard mass distribution: bifocal lens 

In addition to the standard k pdf described above, we also 
try an extreme form for k pdf obta ined using the "bifocal 
lens" mass distribution proposed by iKovn'eil il99 j) : 



p{A) 



S{A-Ai) + 



Ai 



S{A~A2) 



(5) 



Ai- A2 Ai- A2 

where S is the Kronicker's delta function, Ai > 1, and A2 < 
1. For a flxed ko, this pdf can produce more pronounced 
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lensing effects than the standard pdf, if A2 is assigned the 
smallest allowable value, and A\ is made very large. This is 
because the standard pdf has a large probability near A = 1, 
while the bifocal pdf avoids such values altogether. 

The distribution of mass corresponding to eq. |K| is un- 
realistic. One can make it somewhat realistic by allowing 
a range of Ai and A2 values. We do the following: for any 
one line of sight we randomly pick Ai and A2 from specified 
ranges, using flat priors, and eq.|S]sets the probability distri- 
bution of the two amplifications. We use two sets of ranges: 
(1) 1 < ^ (1 -I- 5 Ko) and (1 - kq) < A2 < 1 (we will 
call this model bifocal I model), (2) 1 < Ai s; (1 + 10 acq) 



and (1 



< ^2 ^ 1 (we wiU caU this bifocal II). The 



corresponding n pdfs, obtained by considering many lines of 
sight are shown in Fig. |7[d) as a solid line for bifocal II and 
as a long-dash line for bifocal I. In bifocal II the minimum 
value of A2 is the minimum allowed value; the other limits 
on A's were picked arbitrarily. For comparison, the standard 
K pdf of Section fS. 1.1 1 is also shown, as the short-dash line. 
Optical depth of kq = 0.02 was assumed for all three. Com- 
pared to a standard pdf, bifocal pdfs imply that most of the 
lines of sight are rather empty, and there are a few lines of 
sight with very high values of 



5.2 Biasing scheme 

Projected fractional mass excess, Srimass is related to the ob- 
servationally accessible quantity, Sugai through biasing. We 
chose " power law" biasing, motivat ed by nume rical simula- 
tions of ICen fc Ostriker. a993) and lDekel fcLahav> ClQQSl) . 



(1 + SUgal) = 

[l + Sn 

mass 



1°! +bs (l + <5n„ 
n+b, {l + 5n„ 



for Srlmass > 
for Srlmass ^ 



(6) 



where the power law part allows the biasing to be non-linear, 
with different indexes depending on the sign of Srimass, while 
6s is a stochastic biasing component which is chosen ran- 
domly for each SN; 6s distribution has a Gaussian shape 
and width ai,^. Factor (1 -f Sumass) multiplying 6s ensures 
that the dispersion in Sugai is reduced in underdense regions. 
Qualitatively, the Sugai vs. Srimass relations produced b y 
eq.|S|look similar to those in Fig.l of lDekel fc Lahavl il999tl . 



5.3 Dispersion in the intrinsic SNIa magnitudes 

We assume that the dispersion in magnitude about the per- 
fect standard candle case has a Gaussian shape, with width 

0"SJV- 



6 RESULTS OF MONTE CARLO 
SIMULATIONS 

Each model has five independent parameters: 
{ko, ai, a2, o"i,^ , (TSiv}. We assume flat priors for the three 
biasing parameters, qi, 0:2, and (jf,^, as well as for asN- 
Each specific set of parameters together with one of the 
three mass distribution models generates a synthetic real- 
ization of the 50-source SNIa sample. From this entire en- 
semble of realizations, we only consider those that satisfy 
these observational constraints: (1) the total rms dispersion 
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optical depth probed by galaxies, /c„ 

Figure 6. Probability of obtaining the observed results (slope /3 
of the Mgjv vs. A''> /N relation) given different mass distribution 
models, and a range of optical depths. Solid dots assume standard 
K pdf of Section lS.l.ll empty squares and star symbols represent 
bifocal I and II models, respectively, discussed in Section 15.1.21 
Each point was obtained using 10,000 synthetic realizations of 
SNIa 50-source sample. See Section |H] for details. 



in the synthetic SNIa magnitudes (which includes intrinsic 
and lensing induced contributions) must be within 0.05 mag 
of the actual observed value of O.Smag; (2) the moments 
of the synthetic Sugai distribution (average, standard devi- 
ation, skewness and kurtosis) must be reasonably close to 
those of the observed distribution, which is characterized 
by (0.063, 0.52, 1.7, 4.9). We chose "reasonably close" to 
mean that synthetic values must not be more than a fac- 
tor of 2 away from the actual values. Because of these con- 
straints, the range from which the values of ai, a2, cJb^ , and 
asN parameters are picked is not relevant, as long as it not 
too restrictive. In other words, the observational constraints 
eliminate cases with very large values of these parameters. 

The data suggests that demagnified SNe are not lost 
due to the fiux limit. If they were, we would have less SNe 
at small values of A''> /N, and more at high values of A'^> /N, 
whereas in the data (Fig. the sources are roughly equally 
spread over the entire Ny. /N range. So our synthetic lensing 
models assume that we do not lose SNe because of deampli- 
fication. 

For each synthetic realization of the 50-source SNIa 
sample we derive the corresponding MsN-Ny /N relation. 
The slope of the best-fit line, (3 is recorded. For a given 
value of Ko we generate 10,000 realizations. Solid dots in 
Fig.l^show the percentage of realizations that have (3 smaller 
than the observed value of -0.3, as a function of total opti- 
cal depth, for the standard k pdf. Empty squares and star 
symbols represent the results for the bifocal I and II. 

Figure ^ considers a range of optical depths. What is 
the appropriate value for kq? If the source is at the median 
source redshift for our sample, Zs = 0.47, then APM galaxies 
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Figure 7. Panels (a)-(c) show the dependence of the best-fit 
slope of the Msn^N^ /N relation on the three biasing parame- 
ters (see Section l5.2l . Short-dash lines are for the standard k pdf 
model of Section 15.1.11 while the solid lines are for the bifocal 
II K pdf described in Section 15.1.21 Both models assume optical 
depth Ko = 0.03. The contours are based on 50,000 realizations, 
and are spaced a factor of 3 apart. The shapes of the contours are 
partly determined by the observational constraints, described in 
Section [H] (second paragraph) . It is evident that biasing param- 
eters have little impact on the /3. Panel (d): Short-dash and 
solid lines represent k pdfs corresponding to the contours in pan- 
els (a)-(c) of the same line type. The minimum ft in the bifocal II 
pdf is — Ko (indicated by an arrow). The long-dash line represents 
bifocal I K pdf. 



sample the mass distribution up to a redshift of about 0.3. In 
this case the lensing optical depth probed by APM galaxies 
is 0.011 (left-most arrow in Fig.|SJ. A somewhat more opti- 
mistic estimate for kq is obtained if we assume that galaxies 
trace the mass fluctuations up to the median source redshift, 
0.47, which gives kq = 0.016 (middle arrow). The limiting 
case is obtained by considering the most distant source, at 
Zs = 1.2, and assuming that galaxies probe mass up to a 
redshift of 0.47, which gives kq ~ 0.034 (right-most arrow). 

The important conclusion from this figure is that if the 
standard mass distribution (Section lS.l.ll is assumed, then 
the value of kq does not really matter, and the probability 
of reproducing the observations is < 3%, if kq is within a 
reasonable range. If, on the other hand, the bifocal pdf I or 
II are assumed, the probability of reproducing observations 
depends sensitively on the assumed optical depth. If kq ~ 
0.034 then these models yield 6-20% probability. The bifocal 
pdfs produce more discernible lensing effects, for the same 
Ko because they have a wider range in A, as seen in Fig.EId) . 

The results presented in Fig. |S] marginalize over values 
of all the model parameters (except kq), so the effects of 
the biasing parameters on /? are hidden. However, it turns 
out that these parameters do not correlate with /3. Figure Q 
shows the dependence of /3 on ai, and bs for the bifocal 
II K pdf model (solid contour lines), and the standard k 



pdf model (dashed contour lines), both for kq ~ 0.03. Fifty 
thousand realizations were created for each model, and the 
corresponding contour lines plotted, with adjacent contours 
separated by a factor of 3. It is apparent that for both models 
the dependence on the biasing parameters is weak. 

We conclude that the slope of the MsN-Ny /N relation 
depends on the amount and distribution of matter, i. e. on kq 
and the shape of the k pdf model, but does not depend on the 
specifics of the biasing scheme. Figure |S] quantifies the de- 
pendence on amount and distribution of mass: as expected, 
higher optical depth results in more pronounced lensing ef- 
fects. Standard k pdf produces less lensing for the same kq 
compared to a much broader bifocal pdf. Figure Q demon- 
strates that biasing has little effect on the slope of Msn- 
N-^ /N relation. In other words, the type of lensing signature 
considered here {MsN-Ny /N relation for standard candles) 
probes mass distribution, and is independent of biasing. 



7 CONCLUSIONS 

We detect the signature of weak lensing in the current sam- 
ple of 50 high redshift SNIa taken from two teams: High-z 
Supernova Search Team and Supernova Cosmology Project. 
After corr ecting SNe magnitu des for cosmological distances 
(assuming iTonrv et alj jgjoO?) values), we find that brighter 
SNe are preferentially found behind regions overdense in 
foreground galaxies. This Msn-N^ / N relation has a slope 
P = —0.3 to —0.4, when the angular radii of foreground re- 
gions are 5-15 arcmin. The statistical significance is > 99% 
(see Fig. 01 and 01. The angular radii of 5-15 arcmin imply 
that the lensing structures are 1 — 3 ho.7~^ Mpc across if 
they are located at a redshift of 0.1, and so correspond to 
non-linearly evolved intermediate-scale structure. 

Aside from the possibility that the observed Msn- 
Ny/N relation is a fluke, there is one other possible, non- 
lensing interpretation: Galactic dust obscuration, which 
would make SNIa sources brighter and APM galaxies more 
numerous in the directions devoid of Galactic dust. If this 
interpretation is correct, t hen (1) Galactic ex tinction correc- 
tions applied to the SNe jTonrv et al.ll2003h are much too 
small, and (2) Galactic dust is able to change projected APM 
galaxy number density by 50%, since the observed rms in 
Sugai is ~ 0.5. We cannot comment on (1), but can estimate 
the effect of (2). How much exce ss /deficit in ga l axy d ensity 
can Galactic dust create? Using ISchleeel et all lll99d) data 
we calculate the average Ab Galactic extinction for the 50 
SNe sample to be 0.18 mag. The slope of the R-band APM 
galaxy number counts is dlogN/dm = 0.35, therefore typi- 
cal "dust-induced" {l+Sugai) = lO^B-diosN/dm _ -^g rpj^j^ 
is an upper bound because extinction in the R band will be 
smaller than in ^s, and because extinction averaged over 5- 
15 arcmin radius patches will be smaller than es timates for 
individual SNe, given that ISchleeel et alJ (|^9^ extinction 
map resolution is 6'.1 FWHM. If Galactic dust is indeed the 
cause of the MsN-Ny /N correlation, then Msn and Ny /N 
should separately correlate with SNe extinction. In reality, 
Schle gel ct al.. (.199&) extinction estimates do not show a cor- 
relation with either Msn or Ny/N. Overall, Galactic dust 
interpretation of the observed MsN-Ny /N relation seems 
unlikely, but cannot be ruled out. 

Note that dust intrinsic to the zi ~ 0.1 structures 
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probed by the APM galaxies cannot be invoked to explain 
the observations, because such dust would produce an ef- 
fect opposite to the one detected here, i.e. 13 would be 
positive. If dust is present in groups and clusters traced 
by the APM, it will diminish the amplitude of the effect 
we detect. Presence of dust in groups was suggested by 
iBovle et al.l ||T98S) who found that faint QSO candidates 
are anti-correlated with foreground groups. However, such 
an anti-correlation can also be explained by weak lens- 
ing (iRodrigues- Williams fc Hogai] Il994t ICroom fc Shankj 
Il999r) . Current observations indicate that groups and clus- 
ters do not contain significa nt amounts of dust iFergusonI 
ll993l:lNollenberg et al.li2'ooil . 

If the lensing interpretation in correct, then the stan- 
dard models of mass distribution have some difficulty in re- 
producing P\ the observed value would be detected only in 
< 3% of the cases. We investigate how (3 is affected by the 
amount and distribution of mass along the light of sight to 
the sources, and galaxy biasing schemes. 

We find that larger optical depths and broader k, pdf 
result in steeper /3 slopes (Fig. [SJ. Optical depth is a func- 
tion of global geometry, and no realistic cosmological model 
can give «o ^ 0.05, which is what would be required to com- 
fortably explain the results. Broader k pdfs mean that mass 
fluctuations are more extreme than the standard cosmolog- 
ical models allow, a scenario which is in apparant conflict 
with other means of determining mass fluctuations, like cos- 
mic velocity flows and weak shear lensing. 

We also flnd that biasing has little effect on /3 (Fig. |7|l. 
This insensitivity to biasing is in contrast to weak lensing 
induced QSO-galaxy and GRB-galaxy (anti-)correlations, 
where biasing could, at least partly ex plain the higher 
than expected amplit ude of the effect jjain et alJ l2003t 
I Williams fc Frevll2003l) . We conclude that weak lensing of 
standard candles provides a cleaner probe of the mass distri- 
bution aX z < 0.1, on ~few Mpc scales, then lensing induced 
angular correlations. 

In fact, SNIa can provide the perfect means of mea- 
suring mass inhomogeneities using gravitational lensing. A 
set of standard candles at known redshifts can be analyzed 
using the complementary techniques of weak magnification 
and weak shear lensing. The advantage of magnification lens- 
ing over the more commonly used shear lensing is that with 
the former one can chose the redshift range of the lenses, 
whereas the latter yields the cumulative effect of lensing 
along the entire line of sight to the source. Therefore a large 
uniform set of intermediate redshift SNIa, such as the ones 
that would result from the SNAP mission'^ and the LSST 
mission* would be invaluable for the studies of mass clus- 
tering in the nearby Universe. In addition to detecting in- 
termediate and high-redshift SNe for the purposes of esti- 
mating the global cosmological parameters and the equation 
of state of the dark energy, SNAP will also measure weak 
shear lensing signature duo to large scale structure. Combin- 
ing magnification data (of the type conside red in this paper) 
and shear information llRhodes et alJl2bo3^ for a large set of 
SNe will allow the study of mass distribution at z < 1.5 
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with unprecedented accuracy and 3-dimensional spatial res- 
olution. 

The final issue we address is the impact of weak lens- 
ing on the determination of global cosmological parameters, 
Jim and Ha using SNIa standard candles. In principle, weak 
lensing can affect the derived values of ilm and Qa, if (1) de- 
magnified SNe are preferentially lost from the sample due to 
faint flux cutoff, and/or (2) the k pdf is asymmetric and t he 
SN sample size is small (see also iWambsganss et al.l il997l) ). 
If demagnified SNe were lost from the sample then the dis- 
tribution of SNe in A^> /N would be skewed in the direction 
of larger A'^> /N values. If k pdf is asymmetric, (and it is, 
on scales considered here) , then most SNe in a small sample 
will be slightly demagnified compared to average, and will 
have their Ny /N skewed in the direction of smaller values. 
So both (1) and (2) would make the distribution of SNe in 
A'^> /N uneven, and would to some extent cancel each other. 
In the present sample of 50, the distribution of SNe in A''> /N 
is indistinguishable from uniform, so the average Msn cor- 
responds to the average 7V> /N, and hence lensing effects 
by '^Mpc-size structures probably did not bias the derived 
values of Qm and Qa. 
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